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ABSTRACT 

We examine the effects of dynamical evolution in clusters on planetary systems or 
protoplanetary disks orbiting the components of binary stars. In particular, we look 
for evidence that the companions of host stars of planetary systems or disks could 
have their inclination angles raised from zero to between the threshold angles (39.23° 
and 140.77°) that can induce the Kozai mechanism. We find that up to 20 per cent of 
binary systems have their inclination angles increased to within the threshold range. 
Given that half of all extrasolar planets could be in binary systems, we suggest that 
up to 10 per cent of extrasolar planets could be affected by this mechanism. 
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1 INTRODUCTION 

The majority of star fo rmation is be l ieved to occur in 
clustered environments l|Lada fc Ladal 12001 and refer- 
ences therein). Furthermore, around 60 per cent of solar- 
type stars are observed in binary systems in the field 
ijDuauennov fc Mayor! Il99ll ) . It is believed that a higher 
percentage of these stars (and probably all of them) form 
not as singles, but in bin a ry or higher order sy stems (e.g. 
iGoodwin fc Kroupal [20051 ; IGoodwin et al J 120071 ). although 



Ladal (|2006l ) for an alte rnative interpret ation. 



Se veral studies (e.g. iKroupal Il995al lbl iParker et al.l 
2009a) have demonstrated that intense dynamical evolution 
occurs for stellar systems in typical (i.e. Orion-like) clus- 
ters. If most stars form in such clusters, then this dynamical 
processing may also have a significant effect on any planet 
formation that occurs around binary stars. 

Over 300 extrasolar planets have been discovered to 
dat qj. Of the nearby e xtrasolar planets (i.e. within 200 pc, 
sec iButler et al. 20060, around 40 are known to orbit a 
compon ent of a binary system dDesidera fc Barbierill200"7l ). 
Indeed, iBonavita fc Desideral (|2007 ) suggest that after in- 
completeness has been taken into account, the numbers of 
planets orbiting binary and single stars may be equal. 

Given the significant fraction of extrasolar planets 
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1 The number of discoveries is increasing on an almost daily basis; 
for the most up-to-date catalogue maintained by J. Schneider see 
|http; / / exoplanet.eu/ 1 

^ See http: / /exop lanets.org/planets.shtml for the census by 
Butler et al. (2006) of the extrasolar planets within 200 pc of the 
Earth. 



orbiting the components of binary systems, the dynam- 
ical history of such stellar systems becomes important 
for understanding the pr operties of t he planets therein. 



The Kozai mechanism 
by several authors 



(iKozail Il962h has been 



shown 
also 



e.g. llnnanen et al. I Il997l 
Holman. Touma fc Tremaind Il997l; iTakeda fc Rasid 120051 : 
Malmberg. Davies fc Cha mbers 2007a) to disrupt the orbits 
of planets (see e.g. iMalmberg et~ 2007a, their fig. 1) that 
form in binary systems. 

In this paper we address the possibility of whether the 
dynamical evolution of a cluster may provide a way of induc- 
ing the Kozai mechanism in a significant fraction of binary 
stars. We describe the Kozai mechanism in Section [2] we 
outline our method in Section [3] we present our results in 
Section [4] and we conclude in Section [S] 



2 THE KOZAI MECHANISM 



The Kozai mechanism l|Kozail Il962h was used to quan- 
tify how the orbits of inclined asteroids were influenced by 
Jupiter. It assumes that the mass of the asteroid is negli- 
gible compared to that of Jupiter and the Sun - the same 
assumption can be made for a planet orbiting one of the 
components in a binary system l|lnnanen et al.lll997l ). 

If the inclination angle of the orbit exceeds 39.23°, then 
the Kozai mechanism states that there is a cyclical exchange 
of angular momentum to the asteroid, causing the eccentric- 
ity of the asteroid to vary periodically. The same effect is 
predicted for planetary systems (and protoplanetary disks) 
in binary systems. 

It should be noted that these Kozai cycles can be 
switched off if the inclination angle of the orbit exceeds 
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140.77°. We define a range of inclination angles of systems 
susceptible to the Kozai mechanism as 

39.23° < i Ko z < 140.77°. (1) 

iMalmberg et all (|2007al ) used the Kozai mechanism as 
a method to disrupt planetary systems and showed that it 
induces a highly eccentric orbit for the outer planet, which 
then crosses the orbits of the inner planets. In some cases 
this p rocess leads to t he ejec tion of one or more of the plan- 
ets. |Mahirt)erg^FaI| i|2007ah use intermediate (a ~100AU) 
separation binary systems to invoke this mechanism. 

In our simulations, we investigate whether the inclina- 
tion angle between the components of binary systems can be 
increased from zero (at the birth of the cluster) to an angle 
in the range ikoz (the Kozai angle). We assume that the bi- 
naries at the cluster birth act as gyroscopes and hence any 
systems with an angle greater than the Kozai angle have un- 
dergone significant dynamical processing and could be sub- 
jected to intense perturbations. 



3 METHOD 

3.1 Cluster set-up 

We clo sely follow the method described by IParker et al.l 
(2009a) to set up the clusters and binary systems in our 
simulations. The clusters are designed to mimic a 'typical' 
star cluster, similar to Orion with a mass ~ 10 3 M©. 

For each set of initial conditions, we create a suite of 20 
clusters, identical apart from the random number seed used 
to initialise the simulations. 

We set our clusters up as initially virialised Plummer 
spheres fP lummcr 1911) using the prescription given in 
lAarseth. Henon fc Wielenl l| 19741 ) . The Plummer sphere pro- 
vides the positions and velocities of the centre of mass of 
stellar systems. We adopt three different half-mass ra dii for 
our clusters; 0.1, 0.2 and 0.4 pc. IParker et af] (|2009aT) argue 
that the initial half-mass radius of Orion was in the range 
0.1 - 0.2 pc, but we include simulations with a half-mass 
radius of 0.4 pc for comparison. 



3.2 Binary properties 

We create all our clusters with an initial binary fraction, 
/bin = 1 (i.e. all stars form in binary systems; there are no 
singles or triples, etc.), where 



/bin 



B 



S + B' 



(2) 



and S and B axe the numbers of single and binary systems, 
respectively. 

The mass o f the primary star is chosen randomly from 
a lKroupal (|2002h IMF of the form 



N{M) oc 



M~ 13 m < M/Mq < mi , 
M~ 2 3 mi < M/Mq < m 2 , 



(3) 



where mo = 0.1 M©, mi = 0.5 M , and m2 = 50 Mq. For 
simplicity we do not include brown dwarfs (BDs) in our sim- 
ulations. 

Secondary masses are drawn from a flat mass ratio dis- 
tribution with the constraint that if the companion mass 



is < 0.1 Mq it is reselected, thereby removing the possi- 
bility of choosing secondaries with BD-like masses. This 
constraint maintains the underlying binary fraction, but 
biases the masses of low-m ass systems towards unity (see 
iKouwenhoven et al.ll2009aT l. 

The generating function for or bital periods are the 
log-no rmal distributions obs e rved by iDuquennov fc Mayor! 
1 19911 ') and lFischer fc Marcvl (|l992T i of the form 



/(logP) = Cexp j- 



-(logP-logP) 



(4) 



where logP = 4.8, o"i og p = 2.3 and P is in days. 

Eccentricities of bin ary stars are drawn fr om a thermal 
eccentricity distribution l|Kroupa!ll99"5aL 12008! ) of the form 



Me) = 2e. 



(5) 



Binaries with small periods but large eccentricities 
would expect to undergo the tidal circularisation show n in 
the sample of G-dwarfs in IDuquennov fc Mayor! (|l99ll ). We 
account for this by reselecting the eccentricity if it exceeds 
the following period-dependent value end' 



e tid = ^ (°- 95 + tanh (°- 6P - L7 )) • 



(6) 



This ensures that the eccentricity-period distribu tion 
matches the observations of lDuquennov fc Mayor! (|l99ll ). as 
we expect that tid al circularisation occu rs before cluster evo- 
lution takes place (|Parker et a l. 2009b). Finally, the periods 
are converted to semi-major axes. 

By combining the primary and secondary masses of the 
binaries with their semi-major axes and eccentricities, the 
relative velocity and radial components of the stars in each 
system are determined. These are then placed at the cen- 
tre of mass and centre of velocity for each system in the 
Plummer sphere. 

Simul ations are run using the kira int egrator in Star- 
lab (e.g. iPortegies Zwart et al.l Il999l . l200ll . and references 
therein) and evolved for 10 Myr. 



3.3 Finding susceptible systems 

For each cluster we identify the binary system^] that have 
been preserved since the birth of the cluster (i.e. those that 
have not been broken up by dynamical processing) . We also 
identify the binary systems formed through dynamical inter- 
actions in the clusters as these sy stems are likely to be su s- 
ceptible to the Kozai mechanism (Ma lmberg et al.ll2007bf ). 

Simulations show that disk fragmentation 
is pro bably a major mode of binary formation 



Goodwin. Whitworth fc Ward- Thompson! |2004| ; 



iGoodwin etHI 120071 ) and hence we expect that the orbits 



of the companion and the planets/disk will be roughly 
coplanar (observations of circumstellar disks in binary 
systems have shown that they are usually in clined by less 
than 10 — 20° to the component stars, e.g. Ijensen et alj 



3 We use the neare st-neighbour algorithm described 
by IPa rker et alj ll2009al) (and independently verified by 
IKouwenhoven et alj 12009b ) to determine whether a star is 
in a bound binary system. 
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120041 ; iMonin et al.l 12006), However, this is not an assump- 
tion that affects our resuit as we examine the shift in the 
binary inclination relative to whatever initial inclination 
the planets/disk have. 

We measure the change in the inclination of the binary 
orbit to any planetary system or protoplanetary disk by ex- 
amining the change in the orbital angular momentum vector 
from formation. We assume that the inclination of the plan- 
ets/disk will not be strongly effected by any interaction that 
changes the inclination of the binary and hence any change 
in the inclination of the binary orbit will be a change relative 
to the planets/disk. 

For the preserved systems we use any change in the an- 
gular momentum vector of the binary to ascertain whether 
the system has an inclination angle in the range jkoz after 
each Myr (assuming a birth inclination angle of zero). For 
the systems formed during the cluster's lifetime, we calcu- 
late the change in angular momentum with respect to the 
parental binary or binarieijf] and use the largest angle in our 
calculation^]. 

As we are interested in the effect of cluster evolution 
on planetary systems (or protoplanetary disks), we define a 
periastron distance that allows for the formation of such a 
system without interference from the secondary component 
of the binary system. We take the standard definition of 
periastron distance r pcr i: 

r P cvi = a(l - e), (7) 

where a and e are the semi-major axis and the eccentricity 
of the binary system, respectively. 

We determine r per i for each system, and assume that a 
stable planetary system (or disk) could have developed for 
any binary with r pcr i > 100 AU. Of th e 43 planets orbiting 
a bina ry component in the sample of iDesidera fc Barbieril 
(2007), only 5 are in binary systems with separations less 
than 100 AU. If a binary exceeds this periastron threshold 
we determine its inclination angle from the change in the 
orbital angular momentum vector. 



4 RESULTS 

4.1 The fraction of systems susceptible to the 
Kozai mechanism 

We show the distribution of inclination angles for systems 
with r pcri > 100 AU after lMyr in Fig. □ Whilst only 
1 per cent of all systems have an inclination angle within 
the threshold range (39.23° to 140.77°), many of the binaries 
are tight ( < 30 AU) and hence n ot susceptible to dynamical 
processing l|Parker et alj|2009al ). 

4 We define a parental binary as a system that was a binary 
initi ally and contains one of the stars in a newly formed binary. 
The lDuQuennov &; Mayor] ill 99 if ) period generating function leads 
to very wide (unphysical) binaries in the clusters and so not all 
the primordial binarie s created in the sim ulations are detected 
by our algorithm (see iParker et al.l l2009al . for a more detailed 
discussion). The result of this is that some of the newly formed 
binaries may only have one parental binary. 

5 In practice, in the very few binaries created through these in- 
tense dynamical interactions, both angles tend to be in the range 
!k oz so the angle chosen is irrelevant. 
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Figure 1. The distribution of inclination angles of systems with 
f' P eri > 100 AU after 1 Myr for clusters with an initial half-mass 
radius of 0.1 pc. Binary systems from 20 different realisations of 
the same cluster have been binned together in this histogram. The 
open histogram represents primordial binaries, preserved since 
the birth of the cluster; the hashed histogram represents binary 
systems formed during the dynamical evolution of the cluster. The 
dashed lines indicate the threshold angles (39.23° and 140.77°), 
between which the Kozai mechanism may act upon any planets 
in orbit around the binary stars. 



This is readily demonstrated in Fig. [2] where we plot the 
fraction of binaries with inclination angles in the threshold 
range ikoz as a function of r por i. Binaries with r pe ri > 100 AU 
are far more likely to have inclination angles in the range 
!ko Z than the tighter systems. Additionally, a binary with 
a periastron distance r pcr i < 100 AU is unlikely to form a 
stable planetary system, as that system will be perturbed 
by the other component of the binary system. 

The contribution from binaries that are formed dur- 
ing the cluster's evolution is minimal. As can be seen in 
Fig. [T] only four newly formed systems with r pcr i > 100 AU 
lie within the threshold range «koz, compared to the many 
tens of primordial binaries. 

Of the systems that have r por i > 100 AU, we determine 
the fraction of systems that have an inclination angle be- 
tween 39.23° and 140.77° . We show the evolution of this frac- 
tion during the cluster's lifetime in Fig. [3] We show the frac- 
tion of systems that could be subjected to the Kozai mech- 
anism for three different initial half-mass radii; 0.1 pc (the 
solid line), 0.2 pc (the dashed line) and 0.4 pc (the dashed- 
dot line). 

For the most dense clusters (ri/2 = 0.1 pc), the frac- 
tion of systems that could potentially undergo the Kozai 
mechanism is roughly constant, at ~ 20 per cent. For the 
less dense clusters, the fraction is less; ~ 13 per cent for 
r i/2 = 0.2 pc, and ~ 10 per cent for ri/2 = 0.4 pc. This 
decrease in affected systems with increasing half-mass ra- 
dius is simply due to there being fewer interactions in the 
less-dense clusters. 

The fraction of systems that could be subjected to the 
Kozai mechanism remains roughly constant after 1 Myr for 



4 R. J. Parker and S. P. Goodwin 



0.6 



0.4 



0.2 - 



0.01 




0.1 



1 10 100 1000 

Periastron Distance (All) 



10" 



Figure 2. The fraction of binaries in the range ikoz as a function 
of pcrisatron distance r pcr ; after 1 Myr of cluster evolution (initial 
half-mass radius, rjya = 0.1 pc). 
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Figure 4. The distribution of Kozai timescales, tk oz , in Myr, 
for two hypothetical planets. The open histogram is the distribu- 
tion for a planet with Neptune's mass (5 X 10 _5 Mq) and period 
(165 years), and the hashed histogram is the distribution for a 
planet with Jupiter's mass (9 X 10 — 4 Mq) and period (12 years). 
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Figure 3. The fraction of systems with r pcr i > 100 AU and an 
inclination angle within the Kozai threshold range «Koz (39.23° 
— 140.77°) as a function of time. We show the results for three 
different initial cluster half-mass radii; 0.1 pc, 0.2 pc and 0.4 pc. 
For the most dense cluster, the fraction of systems that could 
be subjected to the Kozai mechanism remains at ~ 20 per cent 
during the first 10 Myr of cluster evolution. 



the most dense clusters (r 1(/2 = 0.1 & 0.2 pc). This is entirely 
due to the crossing t ime of the cluster being only ~ 0.2 Myr 
( Parke r et al . 2009a) and therefore the cluster has relaxed 
and the majority of binaries have already reached dynamical 
equilibrium. 



4.2 The Kozai timescale 

A planet orbiting the component of a binary star 
will undergo Kozai cycles on the following timescale 
( Kis e eva. Eggleton fc Mikkolalfl99i : iTakeda. Kita fc Rasiol 
|200S| ; IVerrier fc Evansll2009h : 



3tt P v 



(1 



2 +m 2 + in p 

1712 



(8) 



where Pbm is the period of the binary, P p is the period of 
the planet, e hin is the eccentricity of the binary, mi and m2 
are the masses of the primary and secondary components of 
the binary respectively, and m p is the mass of the planet. 

For the densest clusters (ri/2 = 0.1 pc), we use Eqn. [8] 
to determine the Kozai timescale for two hypothetical plan- 
ets orbiting a component of each binary. One planet has a 
period and mass similar to that of Neptune (165 years and 
5 x 1O _5 M0 respectively), the other has a period and mass 
similar to Jupiter (12 years and 9 x 10 _4 Mq respectively). 

Our results are shown in Fig. 3] We determine the Kozai 
timescale for the hypothetical planets for all binaries with a 
periastron distance r por i > 100 AU and an inclination angle 
> 39.23°. We show the distribution of Kozai timescales for 
the Neptune analog in the open histogram, and the Kozai 
timescales for the Jupiter analog in the hatched histogram. 

For the Neptune analog, there are a large number of sys- 
tems in which the Kozai timescale is much less than 1 Myr. 
This means that such planets could be subjected to the 
Kozai mechanism very early on in the dynamical evolution 
of the cluster, and given this timescale it is more pertinent 
to ascribe these effects to a protoplanetary disk rather than 
a system of fully formed planets. It should also be noted 
that Kozai cycles could occur several times in 1 Myr for each 
planetary system. 

The Kozai timescale is longer for the Jupiter analog, 
and we would expect on average only one Kozai cycle per 
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Myr for such a planet. However, if a system consists of more 
than one planet, it is the ti mescale for the oute r planet in the 
system that is important (jTakeda et al.l 120081 ). Therefore a 
Jupiter-like planet could undergo more than one Kozai cycle 
per Myr if there were also planets orbiting the star with 
longer periods. 



5 CONCLUSIONS 

We use JV-body simulations to dynamically evolve a typical 
Orion-like star cluster for 10 Myr. We examine the effect 
of dynamical evolution on binary systems that could host 
a stable planetary system or protoplanetary disk. To this 
end we determine the fraction of systems with a periastron 
distance r per i > 100 AU where the secondary component of 
the binary is shifted with respect to the primary by an angle 
within the Kozai threshold range (39.23° < i Ko z < 140.77°). 

For typical clusters with an initial half-mass radius cor- 
responding to that of Orion originally (0.1 - 0.2 pc; see 
iParker et aTll2009al , for a discussion of the primordial half- 
mass radius of Orion), we find that 20 per cent of binary 
systems with r per i > 100 AU have inclination angles in the 
range zkoz that could induce the Kozai mechanism (jKozail 
1 19621 ) . The Kozai mechanism has been shown to drastically 
alter the orbital prope rties of planets l|lnnanen et al.|[l997l ; 
iMalmberg et al.l l2007a|) and in particular the e ccentricities 
of extrasolar planets ( Malmberg fc Daviesll2009l ). 

In these dense clusters, we place two hypothetical plan- 
ets in systems that could be subjected to the Kozai mech- 
anism in order to determine the Kozai timescale, i.e. the 
length of time it takes for a planet to undergo Kozai cycles. 
For Neptune-like planets, we find that the Kozai mechanism 
has the potential to occur more than once a Myr, with a 
longer timescale for Jupiter-like planets. 

For less dense clusters, the fraction of planetary systems 
in > 100 AU binaries that could be subjected to the Kozai 
mechanism is still between 10 and 20 per cent, i.e. a not 
insignificant fraction. 

Around a quarter of the extrasolar planets discovered 
to date are orbiting a compone nt of a binary system wit h 
a separation exceeding 100 AU (|Desidera fe Barbieril 120071 ). 
Due to incompleteness in these observations, it is suggested 
that half of all extrasolar p lanets may be in binary systems 
( Bon avita & De sidera 20Q3) ■ We therefore propose that sim- 
ple dynamical processing of binary stars in clusters could 
feasibly affect up to 10 per cent of all extrasolar planetary 
systems by inducing the Kozai mechanism. 
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